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Least Squares Fit of a Lorentz Function to a set of Experimental 

Points (SADOI Only) 

Complete program 

This routine fits a set of experimental points with the 

function 



f (x) = A + 



B 



i i m 



= A + 



BC 



C 2 + (X-D) 2 




A is the level of the background, B is the height of the 
peak above background, C is half the full width at half the 
height of the peak, and D is that value of X for which f 
attains its peak value. 

Optimum values of the parameters A, B, C, and D are found 
by minimizing the sum of the residuals R: 

R = Z [Y. - f (X.)] 2 W. , 
i=l 111 
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where 

N is the number of experimental points, k < N < 205 :* 
X. is the X-coordinate of the ith experimental point, 
Y. is the Y- coordinate of the ith experimental point, 
W. is the weight associated with the ith point. 

The quantities W. may be: 

1) W. = (l/N) for all i, 

2) W. = (l/N) Z. 

3) W. = (l/N)/Z?~ 

where the Z. are specified by the user. The third alternative 
provides for the weighting of an experimental point by the 
square root of the number of observations or counts at that X, . 

At the start, the user provides the routine with the values 
X. and Y. (and Z. if desired), initial estimates of A, B, C, 
and D, and two parameters, S-, and €. To find the minimum of 
R, the routine proceeds as follows: 

1) In order to achieve greater accuracy, the value of R 
used by the minimization routine is scaled up by some- 
power P of two. Before minimization begins, the routine 
computes R using the given values of A , B-., CL, and D n , 

and tests to see if R. < 2 . If R^ > 2 , the computer 

u — 

stops on an FF002 error stop, which usually indicates 
■ yj 'either that the experimental points cannot be fit by the 

function f , or that the scatter of the experimental points 
:* is so great that the routine may have difficulty in 

determining optimum values of A, B, C, and D. If R Q < 2 

P: is, determihedrsb:i:hat l/l6 < 2? R < l/8,ciahd '' 

p - - ' 0- 
the quantity of 2 I = R r is minimized. 

2) R"' is computed at the points A+S.., B+5.., C + o.,, D + 6, . 
From these values, the approximate direction of maximum 
decrease of R' t is determined. The routine then changes 

A, B, C, and D until a point is reached where the values of 
R "*■ at the points A 1 + S , B + 8 , C + 5 , and D 1 + 5, are 
all greater than the values R' t , at (A, , B , .' C. , D_ )• 
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3) The numbers R 1 , A , B^ C , and D are then printed 
out by the routine. A new value of 8 is computed, 
o 2 = l/lO 5.^ and if S > e, return to step (2) using 
the new 8. 

k) If 8 , < e, the final tolerance, the optimum values of 
R r , A, B, C, and D are^printed, followed by -P as an 
integer. The user can then have the routine compute and 
print f (x) at any desired X, using these optimum values 

of A, B, C, and D. Thus, the final 8, used will be 

-IP 
e < 8 k < lOe. Note that e > 10 . 

Because the same 5 is used for each variable, this routine 

may be deceived into thinking it has found a minimum in 

all four variables. The routine will work best when the 

variation in R r is nearly the same for a change 5 in each 

variable; this may be achieved by scaling the data so that 

A - B 0— 2C o e Also ) it- is required that |A_ + 8-, | < 1, 

f B ± 5 ll < 1 > l C o ± 6 ll < X > and l D ± 5 ll <lm 
The quantity 8 must be chosen with some care. If 8, is 
too large, the non-linearity of the functions f and R 
may cause the routine to hunt in the wrong direction, 
in which case it will eventually hangup; if 8.. is too 
small, the routine will spend too much time searching for 
any minimum at all. 8 should be large enough so that the 
true minimum lies within A Q + 8, , B n + 8 , C n + 8, , 
D Q + 8-, with a fairly high certainty. A goodorule of 
thumb is to take 8 to be the largest of the uncertainties 
in A Q , B Q , G Q , and D Q , or, 1 to 10 percent of the 
value of the largestrof A, B, C, or D. 
DIRECTIONS FOR USE All data is read by routine N12. All numbers are fractions 

< 1, signed, with up to 12 decimal digits. The decimal 
point is assumed to lie immediately after the sign. 
The data tape is prepared as follows; 

1) A single one-hole delay 

2) A title, consisting of any characters except a one-hole delay, 

3) Any number of one-hole delays to terminate the title. 
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The routine will copy all characters between the one-hole 
delays onto the output tape. If no title is desired, 
punch two one-hole delays. 
k) The list of values X., in order , punched for N 12. 

5) A termination symbol, any one of N, J, F, Li 

6) The list of values Y., in exactly the same order as the 

X. . There must be one and only one Y. for each X. . 
1 1 1 

7) A termination symbol, one of N, J, F, L: 

a) If N, W. =. (l/N) for all i. The routine then skips 
to step 10. 

b) If J, F, or L, a list of Z. is to be read, at step 8. 

8) The list of Z.. in the same number and order as the X.. 

1/ 1 

9) A termination symbol, one of N, J, F, L: 

a) If K, W. i= (l/N) Z. . The routine then proceeds to 
step 10. 

b) If J, F, or L, W ± = (i/n)/"z7. 

10) The initial values A , B , C Q , D Q of the variables A, B, C, D. 

11) A termination symbol, one of N, J, F, L: 

a) If N, optimum values will be found for A, B, C, and D. 

b) If J, F, or L, D Q is set to zero, and it is assumed 
that D = in all further calculations . 

12) The values of 5, , the initial mesh size, and €, the final 
tolerance. 

13) A termination symbol, one of N, J, F, or L: 

a) If N, then after the final values of A, B, C, and D 
have been found, the routine will print back X., 
and f(X.) using the final values of A, B, C, and D. 
Then the routine will; proceed to step 1^, 

b) If J, the routine will proceed to step 1^. 

c) If F, the routine will do as in (a) above, then pro- 
ceed to step 16, 

d) If L, the routine skips to step 16. 

Ik) A list of points X.. The computer will print X. and 

J J 

f(X.) as in 13a. The number of points X. may be different 

J J 

from N, the number of points X. . There may be up to 617 

X.'s. 
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15) A termination symbol, any one of N, J, F, or L. 

16) STOP on 2^090. Raising the black switch returns to step 
1, to read in a new data tape. 

STOPS 2^090 raise black switch to read data tape, 

FFOOO Sum check fails on read in 

FF001 wrong number of variables in this list") Skip start 

FF002 Initial sum of residuals exceeds 2~^ J goes to 2^090. 
NOTE In lib it is specified that D should be zero throughout the 

computation. The final results may have D £ 0; this is an 

indication that rounding errors are at least as large as the 

final value of D. 
TIME REQUIRED The running time of this routine is fairly sensitive to the 

initial choices of A, B, C, D, and 5. A rough estimate f or N 

points is ~ 

Time * . N[l +-log, ri — :] seconds. 

_LU <£ 

To this should be added 1 second for each pair of X, f (x) values 
printed out. 
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